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ABSTRACT 

Unitarity of the global evolution is an extremely stringent condition on finite state models 
in discrete spacetime. Quantum cellular automata, in particular, are tightly constrained. 
In previous work we proved a simple No-go Theorem which precludes nontrivial homo- 
geneous evolution for linear quantum cellular automata. Here we carefully define general 
quantum cellular automata in order to investigate the possibility that there be nontrivial 
homogeneous unitary evolution when the local rule is nonlinear. Since the unitary global 
transition amplitudes are constructed from the product of local transition amplitudes, infi- 
nite lattices require different treatment than periodic ones. We prove Unitarity Theorems 
for both cases, expressing the equivalence in 1 -|- 1 dimensions of global unitarity and certain 
sets of constraints on the local rule, and then show that these constraints can be solved 
to give a variety of multiparameter families of nonlinear quantum cellular automata. The 
Unitarity Theorems, together with a Surjectivity Theorem for the infinite case, also imply 
that unitarity is decidable for one dimensional cellular automata. 

KEY WORDS: Quantum cellular automaton; nonlinear dynamics; unitarity. 
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1. Introduction 

Already present in Feynman's inspirational essay on quantum physics and computation [1] 
is the concept of a binary quantum cellular automaton (QCA): a discrete spacetime array 
of quantum processors, each of which has two eigenstates (un/occupied or spin up/down) 
and is coupled to some set of neighboring processors. He explicitly recognizes the diffi- 
culties of reconciling discreteness and locality of interaction with Lorentz invariance — the 
same problem which must be solved in the causal set approach to quantum gravity [2]. 
Grossing, Zeilinger, et al, discovered a similar conffict with group invariance in their at- 
tempts to apply QCAs as quantum simulators, i.e., for quantum rather than deterministic 
computation.* Being unable to reconcile discreteness and locality with translation invari- 
ance and unitarity, they were led to study a class of CAs whose evolution is, although 
'probability' preserving, both nonunitary [6] and nonlocal [7]. 

Near the beginning of our investigation into exactly unitary — and therefore truly 
quantum — CAs we proved that these physically reasonable requirements are actually in- 
compatible for this class of CAs. More precisely, taking the definition of QCA to include 
unitarity, discreteness and locality, we proved [8]: 

NO-GO Theorem. No nontrivial homogeneous linear binary QCA exists on Euclidean 
lattices in any dimension. 

We showed, however, that weakening the homogeneity condition slightly to require 
invariance under only a subgroup of translations allows the existence of linear binary 
QCAs [9]. The simplest of these partitioned [10] QCAs models the evolution of a quantum 
particle and, in the appropriate limit, simulates the 1-1-1 dimensional Dirac equation. This 
physical interpretation motivates two reformulations of the QCA [9] : 

As a quantum lattice gas automaton. The QCA may be formulated as a lattice gas 
with a single particle. Generalizing to multiple particles forces the automaton to be 
nonlinear if we impose an exclusion principle. The dimension of the Hilbert (Fock) 
space of the theory is now exponential in the cardinality of the lattice. If there is no 
particle creation/annihilation, the one particle sector of the Fock space is equivalent 
to the original QCA. 

As a homogeneous linear ternary QCA. In the Dirac equation the amplitudes for 
the particle to be left/right moving at a point are combined into a two component 
field. Equivalently, the QCA may be formulated to have three eigenstates at each cell, 
corresponding to empty, left moving, and right moving. Coupling two copies of the 
original QCA makes the new one homogeneous. 

Both generalizations evade the No-go Theorem; the first by being both nonlinear and not 
quite homogeneous and the second by being ternary. 



* Feynman [3] , Margulis [4] , and more recently, Lent and Tougaw [5] have investigated the possiblilities 
for deterministic computation using QCAs. 
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Emulating Morita and Harao's construction of a computation universal reversible one 
dimensional CA [11], Watrous has recently constructed a one dimensional QCA which 
is quantum computation universal [12] in the sense that it efficiently simulates the uni- 
versal quantum Turing machine defined by Bernstein and Vazirani [13], which in turn 
eflBciently simulates any quantum Turing machine as originally defined by Deutsch [14]. 
Both Morita and Harao's and Watrous' universal CAs are homogeneous, but since each 
may be considered to consist of three coupled partitioned CAs,* they may also be de- 
scribed as partitioned. They are not, therefore, binary — far more eigenstates are required. 
Since a primary motivation for considering QCA models for quantum computation is the 
likelihood that massive parallelism will optimize nanoscale computer architecture [16] and 
the most plausible nanoscale devices have only a few eigenstates [17], it is of interest to 
ask if the No-go Theorem may also be evaded by homogeneous nonlinear binary QCAs. 
In this paper we answer that question affirmatively and, as a first step toward exploring 
the computational power of such architectures, we explicitly parameterize the rule spaces 
for the simplest such QCAs. 

We begin in Section 2 by carefully defining what it means for a CA to be quantum 
mechanical. Rather than restricting a QCA to have only finite configurations as do Watrous 
[12] and Diirr, Thanh and Santha [18], we consider both periodic and infinite lattices. 
In order for the global transition amplitudes to be well defined as the product of local 
transition amplitudes we must distinguish more carefully between these two situations 
than is necessary in the deterministic case. The distinction is emphasized in the following 
section by the observation that infinite QCAs must be asymptotically deterministic, though 
not necessarily quiescent. 

In Section 3 we prove a series of results expressing the equivalence of unitary global 
evolution and sets of constraints on the local transition amplitudes. Utilizing a bijection 
between configurations/pairs of configurations and paths on weighted graphs G1/G2 we 
prove Unitarity Theorems 3.9 and 3.10. These show that unitarity is equivalent to sets of 
constraints on the weights of such paths arising from the condition that the global evolution 
be norm preserving and, in the infinite case, from the additional independent condition 
that the global evolution be surjective; the latter is the content of Surjectivity Theorem 
3.18. 

We observe in Section 4 that the Unitarity and Surjectivity Theorems, together with 
the finiteness of the graphs Gi and G2, prove that unitarity is decidable in one dimension. 
Diirr and Santha have obtained the same result by different methods [19]. Our main 
interest in this section, however, is to extract from these theorems a procedure for finding 
multiparameter families of local rules which define QCAs. We do so, and apply it in the 
cases where the local neighborhood has size 2 or 3. 

At the beginning of Section 5 we show that a pattern observed in the derivation of 



Morita subsequently constructed a simpler universal reversible one dimensional CA which comprises 
two coupled partitioned CAs [15]; presumably the analogous construction also works in the quantum 
context. 
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the small local neighborhood rules in Section 4 generalizes to give families of QCAs for 

any size local neighborhoods. We conclude by discussing connections with deterministic 
reversible CAs and point in directions for further research to explore the computational 
power and physical interpretation of nonlinear QCAs. 

2. Definitions 

A homogeneous CA is defined by a 4-tuple (L, Q, /, E): For the purposes of this paper we 
will take the lattice of cells L to be the integers Z (possibly with periodic identifications to 
Zjv); such a CA is one dimensional. Q is a finite set of states {0,...,^'— 1}; configurations 
are maps a : L ^ Q, the argument of which will be denoted by a subscript. The local rule 
f : Q X C defines the dynamics of the CA which it will be convenient to encode as 

a set of amplitude vectors: for each local configuration A = {ii, ■ ■ ■ ,ik) G Q'^, 

I A)) := (/(0|zi, . . . , ik), ...,/(?- iKi, ... , ik)) e C^. 

We use the variation |.)) of the familiar Dirac notation [20] to denote a vector in while 
distinguishing it from a state vector of some quantum system. The sesquilinear inner 
product on is denoted ((-I-)). Also, the notation for the arguments of / has been chosen 
to evoke that of conditional probability so that /(*| A) is the i^^ component in the amplitude 
vector, given that A is the local configuration. 

If the range of / is {0,1} C C and each amplitude vector has exactly one nonzero 
component, the CA is deterministic: there is a global evolution map F : — > defined 

by 

a'^Fa ^ Vx e L, /(a^la^+s) = 1, (2.1) 

where the local neighborhood E :— {ei,...,eA;} is a finite set of lattice vectors which 
defines the E -sub configurations a^+E '■= (c"x+ei, ■ . . , Cx+Cfe) G Q^. F is well defined since 
for each iJ-sub configuration a^^-^E there is a unique a'^ such that f{o''^W^_^^) — 1. It is 
convenient to take the local neighborhood to be connected; this is no loss of generality 
since the amplitude vectors can be independent of part of the local configuration. In 
one dimension this means the local neighborhood is a sequence of k consecutive integers 
and local configurations can be written as strings ii...ik, ij € Q, indicating states of 
consecutive cells. 

When any of the amplitude vectors has more than a single nonzero component the CA 
is called indeterministic [21]; probabilistic CA models for parallel computation, considered 
already in the original work of von Neumann and Ulam [22], fall into this class. Here 
we are interested in the quantum mechanical situation so the values of / are probability 
amplitudes and the state of the CA is described at each timestep by a configuration vector 
(i){t) E CL*^, the complex vector space with a basis {\a) \ a E L'^} labelled by the 
configurations. The global evolution map F : CL^ —>■ CL*^ is defined on the configuration 
basis by 

F\a) := J2 Fa'al^')^ (2-2a) 
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where 

Fa'a := n fi<K+E), (2.26) 

and then extended linearly to all of the configuration space CL^ . Although the global 
evolution is linear on the configuration space, the local dynamics of the CA is nonlinear 
unless the local rule is constant or additive [23], i.e., unless the amplitude vector is a linear 
function of the local configuration. 

When the CA is deterministic, each factor in the product (2.26) is either or 1. If 
any is the product is 0, otherwise it is 1. Thus definition (2.2) subsumes definition (2.1). 

When the CA is indeterministic, however, we must be more precise about the meaning 
of the product in (2.26). There are two possiblities: 

If L = Zjv there is no problem. All the transition amplitudes F^'^ are defined by 
finite products. To study the CA on an infinite lattice, we may take the usual statistical 
mechanics approach of computing some property [e.g., a correlation function) on a sequence 
of lattices of increasing size and investigating the limiting behaviour of that property. With 
this methodology in mind we refer to a family of CAs (Z^r, Q, /, -E) for all positive integers 
N as a periodic CA and call a configuration admissible for a periodic CA if it is periodic 
with any finite period N. 

If L = Z, the infinite product (2.26) is defined as the limit of a sequence of partial 
products. A nonzero limit only exists if the successive factors in the product converge to 
1. Since the range of / is at most 5^=+^ points in C, this can only happen if beyond some 
stage all the factors in the product are 1. For any configuration cr define the limit set 

n+{a) :={XeQ'' \ yR>0,3x> R: <7^+e = A}, 

and define 0~(a) similarly, i.e., by conditioning on the existence of x < —R. We will 
refer to the two maximal semi-infinite subconfigurations of cr which contain only E- 
subconfigurations from the respective limit set as the ends of a. Then a has a nonzero 
transition amplitude to some configuration only if for each A e Qr^ (cr) some component of 
the amplitude vector |A)) is 1. For any local rule /, let 

Bf — {XeQ^ \ f{i\X) = 1}, 

be the set of local configurations with big amplitude vectors. This description is justified 
by the observation that for each X E Bf, the length of its amplitude vector is greater than 
or equal to 1. Call a set of configurations A C Z*^ admissible iff 

{i) Vcr, a' G A, F^r'a is well defined. 
{ii) "ia e A, 3a' e A such that F^>^ 7^ 0. 
(m) VfT eA,a'eZ^, if F^,^ ^ then a' e A. 

We allow inGnite CAs only when the set of admissible configurations is nonempty. Several 
simple facts follow immediately: 
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Lemma 2.1. A is closed under the evolution. Furthermore, if a E A then Q^(a) C Bf, 
and if n^{a') = n^{a) then a' e A. 

The most familiar infinite CAs have a unique quiescent state [22], say 0, defined by the 
property that /(^|0 . . .0) = 5io, i.e., Bf = {0 . . .0}. (These are the only QCAs allowed in 
[12] and [18,19].) For such CAs admissible configurations have finite support in the sense 
that the only nonquiescent cells lie in some finite domain. It is clear that the collection 
of configurations with finite support is closed under the global evolution and, depending 
on /, may satisfy the admissibility conditions in Section 3 we show that these 

properties can hold for more interesting sets of configurations as well. 

Consider the subspace of CL'^ with basis the admissible configurations A. This will 
become a Hilbert space, the physical configuration space H, once we define an inner 
product. Again using Dirac notation [20], but now for vectors [V') ^ denoting 
states of a quantum system, 

where |0) = Xlo-eA '^'^1'^) similarly for l-i/^). H C CA is the set of vectors with finite 
norm. A periodic or infinite CA is a quantum CA if the global evolution preserves this 
inner product and is time reversible (which is a distinct condition for infinite dimensional 
Hilbert space [24]), i.e., is unitary. The physical interpretation is that when the system is 
described by the configuration vector |0) with (</)|</)) = 1, the probability of observing it to 
be in configuration cr is (p^cj)^. The invariance of the inner product implies conservation of 
probability [20,24]. 

3. The unitarity constraints 

The condition that the global evolution of a QCA be unitary places strict constraints on 
the local rule /. These constraints can be expressed directly in terms of the amplitude 
vectors: 

Theorem 3.1. F is unitary iff for all admissible configurations a' and a", 

n H'^x+eWx+e)) = S<7"a', (3.1) 

xeL 

and, in the infinite case, F is surjective. 

Proof F is unitary iff it preserves inner products and is surjective. For finite dimensional 
F the former condition implies the latter, as well as that F~^ = F^. In the infinite case 
both conditions are necessary and sufficient to reach this conclusion [25]. The condition 
that F preserves inner products is equivalent to / = F^F; in components this becomes: 

Sa"a' — ^ (-^^)<t"ct-PV<7' = ^ Faa"Faa' 
a£A aeA 
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= En/KK+^)n/KK+^) 

a€Ax£L y£L 

= E n fKK+E)fM^UE) 

(where the eUipses denote sums for each Xi & L and by admissibihty condition {Hi), only 
a & A make nonzero contributions) 

= n E fKK+E)fM<rUE) 

xEL cr^eQ 

= YI{{<^x+e\<^'x+e)), 
x£L 

which is, of course, no more than the famihar condition for unitarity we described in 
[8] from the perspective of the sum-over- histories approach to quantum mechanics [26], 
written in terms of the transition amplitudes (2.26). ■ 

A priori Theorem 3.1 tells us only that inner product invariance requires a constraint 
for each pair of admissible configurations. To be able to verify unitarity for arbitrarily large 
or infinite lattices given a local rule we must reduce the number to something independent 
of the cardinality of L. We do so in Subsections 3.1 and 3.2; then in Subsection 3.3 we 
consider the surjectivity condition which must also be satisfied for infinite lattices. Our 
results in this section are similar to those obtained independently by Diirr, Thanh and 
Santha [18,19]. 

3.1. The normalization constraints . . . 

There is one diagonal constraint in (3.1) for each admissible configuration: 

l=\{{{a^+EW.+E)). (3.2) 

x£L 

That is, for each admissible configuration, the product of the lengths of the amplitude 
vectors of its £'-subconfigurations must be 1. So to understand the structure of this set 
of constraints, we must understand which sets of £^-subconfigurations can occur together. 
Since our configurations are one dimensional they may be constructed by appending (and 
prepending) one element of Q after another; each successive element is the right (left) most 
component of a corresponding £'-subconfiguration which is determined by the new element 
together with the previous k — 1 elements. 

This construction may be realized by a directed graph (known as a de Bruijn graph 
[27]) which is weighted [28]: Let Gi{Q, k) be the (g, g)-valent directed graph with vertices 
labelled by the elements of Q^~^ and directed edges labelled by Zfe e Q connecting vertices 
labelled ziZ2 . . - ik-i to vertices labelled Z2 . . .ik-iik- (See Figures 1 and 3 in Section 4 for 
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examples.) A path of length / in Gi((5, k) is a sequence of I directed edges in the graph, 
each of which (but the first) is directed away from the vertex to which the previous edge is 
directed. Paths of length / correspond bijectively to subconfigurations of length I + k — 1; 
each edge corresponds to an £'-subconfiguration.* A cycle in a directed graph is a closed 
path which passes through no vertex more than once. 

Assign the weight {{iii2 ■ ■ -ikliih ■ ■ -ik)) to each edge labelled ik leaving a vertex 
labelled iii2 ■ ■ ■ ik-i and define the weight of a path in the directed graph to be the product 
of the weights of its edges. Then the weight of the path corresponding to configuration a 
is given by the right hand side of (3.2). 

Lemma 3.2. If (3.2) holds for every admissible configuration then the weight of every cycle 
in G\{Q, k) is 1. 

Proof. Since any cycle corresponds to an admissible configuration cr of a periodic CA for 
some A^, (3.2) must hold for a, proving the statement in the periodic case. In the infinite 
case, consider any cycle with edges labelled and a bi-infinite path in Gi{Q,k) 

containing the cycle which corresponds to an admissible configuration cr. cr must have a 
subconfiguration corresponding to the cycle starting after some Xq G Z: <7xo+j = ij for 
1 < j < I and a^o+i+j = ij for < j < k — 1. Define a new configuration a' by 

^ { cr^ if X < xo; 
* \ a^+i ii x> Xq. 

O^(cr') = O^(cr), SO a' is admissible by Lemma 2.1. Furthermore, the product in (3.2) will 
be identical for a and a' but for an extra product in the former of the weights of edges in 
the cycle. Since (3.2) holds for both cr and cr', the weight of the cycle must be 1. ■ 

Theorem 3.3. For a periodic CA, (3.2) holds for all admissible configurations iff the 
weight of every cycle in Gi{Q,k) is 1. 

Proof. Lemma 3.2 is the 'only if part of this statement. The 'if part is almost a tautology 
since the only configurations in a periodic CA correspond to collections of cycles, so if each 
has weight 1, then (3.2) holds for all admissible configurations. ■ 

The analogous result in the infinite case is somewhat more complicated because the 
admissible configurations are different. Let us examine them more closely: 

Lemma 3.4. Let a be any configuration in an in&nite CA. Then Qr^{a) consists of E- 
subconhgurations corresponding to the edges in a collection of cycles in Gi{Q, k). 



* Equivalently, Gi{Q, k) is a finite state automaton (FSA; see, e.g., [29]) with states corresponding to 

length k — 1 subconfigurations and transitions to single elements of Q. Such FSAs model circuits 
with memory, like convolutional coders [30]. When the outputs of an FSA are (state, transition) pairs 
it is called a Mealy machine [31]; here the outputs are -B-subconfigurations. 
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Proof. We must show that for any £'-subconfiguration A e fl'^{a) there is some cycle in 
Gi{Q, k) containing the edge corresponding to A such that the i?-subconfiguration corre- 
sponding to each edge in the cycle is also in Q~^{a). By definition, A G Q~^{a) only if the 
corresponding edge in Gi{Q, k) is traversed infinitely many times by the infinite forward 
path corresponding to cr. To return to this edge the infinite forward path must contain a 
cycle; to return to it infinitely many times it must contain at least one of the only finitely 
many cycles infinitely many times. Each edge in that cycle is therefore traversed infinitely 
many times by the infinite forward path corresponding to a and hence the correspond- 
ing iJ-subconfigurations are in Vf^^a). The analogous proof with 'forward' replaced by 
'backward' proves the statement for f2~(cr). ■ 

Applying Lemma 2.1, we have the immediate: 

Corollary 3.5. An inGnite CA has admissible conGgurations only if there is a subset 
Df C Bf of local configurations corresponding to the edges in a collection of cycles of 
Gi{Q,k). 

Let Di{Q,k,f) be the subgraph of Gi{Q,k) consisting of these cycles. A maximal 
Df which is closed under the evolution will be referred to as the deterministic sector of a 
QCA because of the following result: 

Theorem 3.6. If F is unitary then when restricted to local configurations in Df, the local 
rule is deterministic. If Df is closed under the evolution then the set of configurations 
defined by A := {a \ ^l^{a) C Df} is admissible. 

Proof. By Corollary 3.5, each X & Df G Bf corresponds to an edge in a k) cycle, 

each edge of which corresponds to a local configuration Xi E Df, 1 < i < I. By Theorem 
3.1 and Lemma 3.2, if F is unitary. 

When we defined Bf we observed that for each A G Bf, |A)) has norm at least 1; now 
(3.3) implies that each Aj in this cycle, and hence every local configuration in Df, has 
an amplitude vector with norm exactly 1. Since each amplitude vector has at least one 
component equal to 1, all the remaining components must be 0; thus the local rule is 
deterministic. 

Since F restricted to local configurations in Df is deterministic, for any a, a' G A, 
either the ends of a evolve to those of a', in which case F^'cr is well defined; or they do not, 
in which case -Fo-'cr = 0, which is still well defined. This verifies admissibility condition (i). 
Furthermore, Df is closed under the evolution means that for every path of length k in 
Df, corresponding to a subconfiguration ii . . . i2k-ii the unique local configuration i'^ . . . 
defined by /(^^Kj • • •■ij+fc-i) = 1, 1 < J < ^, corresponds to an edge in Df. If F is unitary 
then it is norm preserving, so for any a E A, there must be some a' G with Fo-'o- ^ 0. 
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Since Df is closed under the evolution, Cl^{a') C Df, this verifies admissibility conditions 
(ii) and, since Df is maximal, [iii). ■ 

Since we only consider infinite CAs with admissible configurations, we will assume 
henceforth that each has a nonempty deterministic sector which defines the set of admissi- 
ble configurations as in Theorem 3.6. The remarkable consequence of this theorem is that 
an infinite QCA must be asymptotically deterministic: the only admissible configurations 
are those which evolve deterministically outside some finite domain. Indeterministic quan- 
tum evolution can only occur for a finite subconfiguration interpolating between the ends 
of the configuration. Thus the analogue of Theorem 3.3 for infinite CAs is: 

Theorem 3.7. For an inGnite CA, (3.2) holds for all admissible conGgurations iff the 
weight of every cycle in Gi{Q,k) is 1 and the weight of every acyclic path in Gi{Q,k) 
terminating at vertices in Di{Q, k, /) is 1. 

Proof. By Lemma 3.2, if (3.2) holds for all admissible configurations then the weight of 
every cycle in Gi{Q,k) is 1. Consider any acyclic path connecting cycles in Di{Q, k, f). 
By Theorem 3.6, there is an admissible configuration corresponding to a bi- infinite path 
contained in -Di(Q, k, f) except for a finite segment along the chosen acyclic path. If (3.2) 
holds for all admissible configurations the weight of the whole path is 1, as are the weights 
of the two ends in Di{Q, k, /); hence the weight of the acyclic path is also 1. Conversely, 
by definition and by Lemma 3.4, the only admissible configurations are those with ends 
corresponding to collections of cycles, connected by a finite subconfiguration which, after 
removing some finite number of subconfigurations corresponding also to cycles in Gi{Q, k), 
is an acyclic path terminating at vertices in Di{Q, k, f). If all the cycles and the acyclic 
path have weight 1 then (3.2) holds. ■ 

3.2. ... and the orthogonality constraints 

There is one off-diagonal constraint in (3.1) for each pair of distinct admissible con- 
figurations: 

0=11((<+eK+e))- (3-4) 

By (3.2), for each x such that o-'J._^_^ = o-'^^^, the corresponding factor in (3.4) is nonzero. 
Thus at least one of the pairs of mismatched £'-subconfigurations o-'J._^^ ^ ^^'x+e ni^st 
contribute a factor of in order for the product in (3.4) to vanish. To understand this set of 
orthogonality constraints, therefore, we must understand which sets of pairs of mismatched 
^^-subconfigurations can occur together. We construct a new weighted directed graph 
which generates these sets: 

Let G2{Q,k) be the (g^, g^) -directed graph with vertices labelled by the elements of 
Q^~^ X Q^~^ and directed edges labelled by (i'l^i'i^) G connecting vertices labelled 

{iii2 ■ ■ ■ i'k-n'i-'A ■ ■ ■ ■^fe-i) to vertices labelled (^2 • • • 'i-'k-i^k-i - ■ ■ '^'k-A)- Note that there 
is a subgraph of G2{Q,k) isomorphic to Gi{Q,k), namely those vertices and edges with 
both components of these labels identical; so we write Gi{Q, k) C G2{Q, k). (See Figure 2 
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in Section 4 for an example.) 

Paths in G2{Q, k) correspond bijectively to pairs of subconfigurations; each edge cor- 
responds to a pair of £'-sub configurations. We assign the weight ((^l^2 ■ ■ - "^kVA - ■ ■ ^'kfl 
each edge labeUed leaving a vertex labelled (^l^2 • • ■i'l-ii'^'A ■ ■ -^fc-i) ^^^^ ^^e 

weight of the path corresponding to a pair of configurations a' and a" is given by the right 
hand side of (3.4). Two configurations a' and a" are distinct only if they have at least one 
pair of mismatched E'-subconfigurations, so a path in G2{Q, k) corresponds to them only 
if it intersects M{Q,k) := G2iQ,k)\Gi{Q,k). 

Lemma 3.8. If F is unitary then the weight of every acycUc path in M{Q, k) terminating 
at vertices of Gi{Q, k) is 0. 

Proof. Any acyclic path in M{Q,k) terminating at vertices of Gi{Q,k) can be extended, 
using only edges in Gi{Q, k), to a path in (?2(<5, k) corresponding to a pair of distinct ad- 
missible configurations a' and a". If F is unitary, however. Theorem 3.1 implies condition 
(3.2), which precludes the weight of any edge in Gi{Q, k) from vanishing. Thus (3.4) can 
only be satisfied for a' and a" if the weight of the portion of the path in M (Q, k) vanishes. 

■ 

Theorem 3.9. For a periodic CA, F is unitary iff all of the following are true: 

(i) The weight of every cycle in Gi{Q, k) is 1. 

(ii) The weight of every cycle in M{Q,k) is 0. 

(in) The weight of every acyclic path in M{Q, k) terminating at vertices in Gi{Q, k) is 0. 

Proof. If F is unitary then Theorem 3.1 and Lemma 3.2 imply (z) while Theorem 3.1 and 
Lemma 3.8 imply (Hi). Any cycle in M{Q, k) corresponds to a pair of distinct admissible 
configurations a' and a" of a periodic CA for some N. Theorem 3.1 implies (3.4) must 
hold for a' and a"; this implies {ii). 

Conversely, by Theorem 3.3, {i) implies (3.2) holds for all admissible configurations in 
a periodic CA. Furthermore, pairs of distinct admissible configurations in a periodic CA 
correspond to collections of cycles in G2{Q,k), at least one of which intersects M{Q,k). 
Thus (3.4) holds if each cycle which intersects M{Q, k) has weight 0; [ii) and (Hi) imply 
that this is the case. By Theorem 3.1, (3.2) and (3.4) imply F is unitary. ■ 

Before stating the analogous result in the infinite case, we define D2{Q,k, f) to be the 
subgraph of G2{Q, k) consisting of the edges corresponding to pairs of local configurations 
each of which is in Df. 

Theorem 3.10. For an inRnite CA, F is unitary iff all of the following are true: 

(i) The weight of every cycle in Gi{Q, k) is 1. 

(ii) The weight of every acyclic path in Gi{Q, k) terminating at vertices in Di{Q, k, f) is 
1. 

(Hi) The weight of every acyclic path in M{Q, k) terminating at vertices in Gi{Q, k) is 0. 
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(iv) The weight of every cycle in D2{Q, k, f) D M{Q, k) is 0. 

(v) F is surjective. 

Proof. If F is unitary then Theorem 3.1 and Theorem 3.7 imply (i), (m) and (v), while 
Theorem 3.1 and Lemma 3.8 imply (m). Any cycle in D2{Q.,k,f) fl M{Q,k) can be 
repeated infinitely and then corresponds to a pair of distinct admissible configurations a' 
and a" (since the ends of each are contained in Df). Theorem 3.1 implies (3.4) must hold 
for a' and a"; this implies {iv). 

Conversely, by Theorem 3.7, (i) and (n) imply (3.2) holds for all admissible config- 
urations in an infinite CA. Furthermore, pairs of distinct admissible configurations in an 
infinite CA correspond to bi-infinite paths in G2{Q, k) which consist of cycles in D2{Q, k, /) 
at both ends. If any one of these cycles lies in D2{Q,k,f ) fl M{Q,k) then {iv) implies 
(3.4) holds. If none of the cycles contains only edges corresponding to mismatched E- 
sub configurations then there must be an acyclic path in M{Q, k) terminating at vertices 
in Gi{Q,k). In this case (m) implies (3.4) holds. By Theorem 3.1, (u), (3.2) and (3.4) 
imply F is unitary. ■ 

3.3. Surjectivity 

To complete the program of determining a collection of constraints on the amplitude 
vectors which is equivalent to unitarity of the global evolution we must show that for 
infinite CAs, condition {v) of Theorem 3.10 — surjectivity — is equivalent to some collection 
of such constraints. 

By Theorem 3.6 each a E A has a minimal length interior subconfiguration I{a) = 
(7x0+1 ■ • • o'xo+n, defined by the condition that if j < xq + 1 or j >xo + n + k — 1 then 
(Tj-k+i . . . CTj e fl^{(7), respectively. We will refer to the semi-infinite subconfigurations 
. . .a^o and axo+n+i • • • as the deterministic ends of a. For example, if Df = {000, 111} 
{k = 3) then a possible admissible configuration is a = . . . Oai . . . a^l . . . and I{a) = 
(Ti . . .(7„, provided cri ^ and (j„ 1. Let := {a e A \ |/(cr)| = n}. Then A = UA„, 
< n < oo. We will find a necessary and sufficient set of constraints on / such that F 
maps onto A^ for each nonnegative integer n; since A is an orthonormal basis for H this 
is equivalent to F mapping onto H. 

We begin by observing that F already maps onto the completely deterministic con- 
figurations Aq. More precisely. 

Lemma 3.11. For an infinite CA, if F is norm preserving (i.e., if conditions (i) — (iv) of 
Theorem 3.10 hold) then F : Aq ^ Aq is surjective. 

Proof. If F is norm preserving it is injective on H D Aq. By definition, all i?-sub configura- 
tions of Aq are in Df, so F is deterministic on Aq.* Every configuration a E Aq corresponds 



* Notice that neither Moore and Myhill's Garden of Eden Theorem [32] nor Hedlund's result that 
injectivity implies surjectivity for endomorphisms of shift dynamical systems [33] may be applied to 
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to a sequence of cycles in Di{Q, k, f) since I{a) = 0. But injectivity implies surjectivity 
for F restricted to cycles in Di{Q, k, /): Each cycle of length m must map to a closed path 
of length m which, by injectivity, must also be a cycle. There are only a finite number of 
length m cycles and by injectivity F maps no two to the same one; hence F is surjective 
on cycles in Di{Q, k, f). ■ 

This result has an immediate corollary once we make some suitable definitions: 

Corollary 3.12. For an infinite CA, if F is norm preserving then the restriction of F to 
semi-infinite completely deterministic admissible subconGgurations is surjective. 

To define the restriction of F on subconfigurations it will be convenient in the rest of 
this subsection to take neighborhoods of size k to be defined hy E := {0, . . . ,k — 1}. Any 
other connected local neighborhood of the same size is equivalent to this choice after a 
translation, so there is no loss of generality Having made this choice for E, we may define 
the restriction of F to any subconfiguration ... cr^ of any admissible configuration by 

b-k+l 

F\aa...(Tb) ■■= K---<^b-k+i) n /(^*Vi---f^i+fc-i)> (3-5) 



which is well defined even when either a = — oo or 6 = oo since a"a . . . a"b is a subconfigura- 
tion of an admissible configuration. ((3.5) is equivalent to (2.2a, b) when a = — oo, b = oo.) 
Now we can prove Corollary 3.12: 

Proof. By definition, a semi-infinite completely deterministic admissible subconfiguration 
is a deterministic end of at least one configuration in ^o- By Lemma 3.11, F is surjective 
on ^o; whence the restriction of F defined by (3.5) to subconfigurations of configurations 
in ^0 is also surjective. ■ 

This means that F restricted to the deterministic ends of admissible configurations is 
surjective and indicates that we may concentrate just on interior subconfigurations. For 
A = Ao . . . Afc-i, p = po • • • Pk-i e Df define 

"^n^ '■= {n ■ ■ • inPO ■ ■ ■ Pk-2 I «i e Q, 1 < j < n} 

^(A,p) ._ ^x^_ Xk-iH ■ ■ ■ inPo ■ ■ ■ Pk-2 \ ij ^ QA < j < n}. 

Considering an element of Ai^'^^ to be a subconfiguration of an admissible configuration 
we may use (3.5) to define the action of F on A^'''\ 

Theorem 3.13. For an infinite CA with F norm preserving, F maps onto An for < n G Z 
iff for all X,p' e Df, its restriction F^^'^ -* : A^'^^ —>■ An is surjective when f{p'^_i\p) = 1. 



conclude that F is surjective on Aq since A, and hence Aq, is only a proper subset of the set of all 
possible configurations. 
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Proof. Each configuration a' G has an interior subconfiguration of length n extended by 
deterministic ends: a' — . . . X'i^_ii'i . . . i'nP'o ■ ■ ■ where a'^^^ = i'l, say, and X'j_j^_^_i • • • A^- & Df 
for j < k — 1, Pj . . . Pj_^_i._i e Df for j > 0. F maps onto iff for all a' e A^, 

\a')^Fj2cmWn (3.6) 

for some element Cm\o''^) € H. By Corollary 3.12, there exist semi- infinite deterministic 
configurations . . . Afc_i and po . . . such that {a'\F\a^) ^ only for configurations of the 
form = . . . Aq . . . Xk-i^T ■ ■ ■ CPo • • •> where a^^ = and 



F| . . . Aq . . . X,^_^) — I . . . X'i^_^) 
F\po---) = Ipq...). 



(3.7) 



For these configurations (3.5) implies that 

= {i[ . . . i'J^ . . . p'k_2\F\X, . . . X,_,iT . . . CPo • • • Pk-2)- (3.8) 
Combining (3.6), (3.7) and (3.8) we find 



= Fj2cm\X,... X,_,iT . . . CPo • • • Pk-2), (3.9) 



which is exactly the statement that Fn is surjective. Since A' := Ao...A^_j^ G Df 
is arbitrary and /(A'^_]^|A) = 1 (from (3.7)) defines a unique A := Ao...Afc_i G Df by 
Corollary 3.12, A is an arbitrary deterministic local configuration. Furthermore, (3.7) 
implies f{p'j^_i\p) = 1 for p := po---Pfe-i; this completes the "only if" direction of the 
proof. Conversely, if for all X, p' G Df with f{p'i^_i\p) = 1, (3.9) holds for some Cm G C, 
then using Corollary 3.12 we can construct arbitrary deterministic ends satisfying (3.7) 
and conclude that (3.6) holds for all a' G A^- ■ 

With Lemma 3.11 and Theorem 3.13 we have reduced the problem of the surjectivity 
of the infinite dimensional map F to the equivalent problem of the surjectivity of the 
finite dimensional maps Fn for all 1 < n G Z (we will supress the superscripts (A, p') in 
the following). But F^ is onto iff detF^ ^ 0. To understand these conditions we must 
explicate the structure of F„. Notice first that as an immediate consequence of definition 
(3.5) each column of Fn has a common factor. More precisely. 

Lemma 3.14. For a, a' G Q"^, 

(P^ / (P0---Pfc-2l^l«n-fc+2---«nP0---Pfe-2): if U > k - 1; 

{^n)a',a - {^n)a',a " j (p^ . . . p',_,\F\Xn+, . . . X,_,ap, . . . p,_,), if < Tl < k - 1; 

(3.10) 

where F is the reduced transition matrix defined by 

iFn)a',a := {a'\F\Xi . . . Xk-ia). (3.11) 
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These common factors can be pulled out of each column in the determinant: 
Corollary 3.15. detF^ = Cn detF^, where 

/n-y6Qfe-i(Po---Pfc-2l^l7Po---Pfc-2)'^" ifn>k~l; 
YlaeQ-(Po---Pk-2\F\K+i---K-iO^Po---Pk-2)^ ifO<n<k-l. 



Proof. When n > A; — 1, for each 7 e Q^~^ there are columns in F„ labelled by- 

sub configurations having 7 as their last k — 1 states. According to Lemma 3.14, the factor 
in (3.10) occurs in each. When < n < /c — 1, for each a G there is exactly one column 
in Fn labelled by a. Again, by Lemma 3.14, this column contains the factor in (3.10). ■ 

Second, for 7 e Q'^~^, let 

/ I7O)) \ 



$(7) 



e Mq{C). (3.13) 



\h{q-l)))J 



Now we can find a recurrence relation (with initial condition Fq — 1) for the reduced 
transition matrix: 

Lemma 3.16. -Fn+i can be expressed in terms of Fn and the ^^'^^ as 

where a, a' G and i^j G Q. 
Proof Apply definition (3.11): 

{Fn+l)a'j,ai = (a'j|-^|Ai . . . Xk-iai) 

= {a'\F\X,...X,_,a).i{[f--'+'\-''-'\ ^! " 7 ^ " !' , 
^ ' ' ^ \/(j|An+i...Afc_im), ifO<n</c-l; 

where the second equality follows from definition (3.5). Definition (3.13) gives the result 
(3.14). ■ 

Just as the common factors in Lemma 3.14 could be pulled out of the determinant of 
Fn in Corollary 3.15, so too can the determininants of the common tensor factors 
We need two simple properties of determinants: 
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1. For i,j e {l,...,m}, let e Mn{C) he m x m matrices over C. Suppose 
Xa = XiB with 0;^ e C and 5 e M„(C). Then 

(Xil X12 ■ ■ ■ Xim \ 

: : : h (3-15) 

for / the n X n identity matrix. 

2. Let X e M^(C). Then 

det(X /) = (det X)" (3.16) 
where / is again the n x n identity matrix. 

Corollary 3.17. detF„+i = dn+i[det Fn]'^ , where 

:H n^^«'-'''^r''r"", (3.17) 

^ naeQ"det$(^"+i -^^-i'*), ifO<n<A;-l. ^ ^ 



Proof. By Lemma 3.16 the columns of come from the columns of tensored with 
the corresponding Applying (3.15) we have: 

det F„,. = det F„ « / . I n.e«'- . [det if n > . - 1 ; 

^ \naeQ"det$(^"+i-^^-i«), ifO<n</c-l 

(where / is the q x q identity matrix), by exactly the same reasoning as in the proof of 
Corollary 3.15, but now using the result of Lemma 3.16 for the common tensored matrices. 
Applying (3.16) gives the result (3.17). ■ 

Putting these results together gives the Surjectivity Theorem summarizing the condi- 
tions that for all 1 < n e Z, det F„ ^ 0: 

Theorem 3.18. An infinite CA with norm preserving F and deterministic sector Df is 
surjective iff for all A, p, p' G Df such that f{p'j^_i\p) = 1, none of the following vanish: 

(i) {p'o...p'k_2\^hpo---Pk-2) 

(a) det$(T) 
for any 7 G Q^~^ , < n < A; — 1. 

Proof. By Corollaries 3.15 and 3.17, detF^ 7^ for all n > iff none of the defined in 
(3.12) nor the dn defined in (3.17) vanish. This is ensured by (i) and (ii), respectively, 
since 

{\Xn+i . . . Xk-iapo . . . Pk-2) \ aeQ'',0<n<k-l} G {\^po . . . pk-2) \ 7 e Q^~^} 

and 
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|^(A„+i...Afe_ia) I CK e Q^, < n < /c - 1} C {$(^) I 7 e Q''-^}. 



We will refer to the nonvanishing of the expressions in (i) and (m) as the surjectivity 
constraints. Each is a closed condition, so as we will see in the next section, when they 
are not inconsistent with constraints {i)-{iv) of Theorem 3.10, they do not reduce the 
dimension of the solution space. 

4. Solutions 

Theorems 3.9 and 3.10, to which we will refer as the Unitarity Theorems, and Surjectivity 
Theorem 3.18, together with the finiteness of the graph G2{Q, k), show that given a local 
rule / we can determine, by checking only a finite number of conditions, whether the global 
evolution is unitary, i.e., whether / defines a QCA. This proves: 

Theorem 4.1. Unitarity is decidable for 1 dimensional CAs. 

As yet we have no proposed local rule / for which to check unitarity. In the next two 
subsections we write down the constraints resulting from the Unitarity Theorems for binary 
{i.e., q = 2) CAs in the two simplest cases and show that they can be solved to give 
multiparameter families of QCAs. 

Before doing so, however, it is useful to discuss some symmetries of QCAs. First, note 
that the symmetry group of a one dimensional lattice L has two generators: T, translation 
by one and P, refiection in the origin/parity reversal. P acts on local configurations, 
sending A = zi . . . to PA = ik ■ ■ - ii- (To be precise, this is the case when k is odd; when 
k is even this transformation is TP if the origin is taken to be at the |_/c/2j position in the 
local neighborhood. It will cause no confusion, however, to denote both by P.) P may be 
defined to act on a local rule / by: {Pf){i\X) :— f{i\PX); hence P\ii . . . ik)) := \ik ■ ■ - h))- 
We will refer to Pf as the parity transform of /. Since P^f = f, any local rule which is not 
symmetric, i.e., invariant under this reflection, will pair with a distinct parity transformed 
local rule. 

Second, note that the symmetric group Sg acts on the set of states Q of a QCA. For 
q = 2 the symmetric group is generated by the transposition r which interchanges and 
1.* T acts on a local conflguration \ — ii . . .i^ to give tA = rii . . . rik, so it acts on a local 
rule in two ways: 

(Ti„/)(z|A) :=/(z|tA) 
(Tout/)(^|A) :=/(tz|A). 

These transformations will be most useful in the discussion of inflnite, and hence asymp- 
totically deterministic, QCAs. 



* Prom a more physical perspective, r is analogous to charge conjugation. 
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4.1. k = 2 

The simplest nontrivial binary CAs 
have local neighborhoods of size 2* Fig- 
ure 1 shows Gi. (The arguments Q = 
{0, 1} and k = 2 are suppressed in this 
subsection.) Here all the relevant paths 
can be identified by inspection. More for- 
mally, define the transfer matrix Ai for 
Gi to be the 2x2 matrix with ij^^ en- 
try the weight of the edge from vertex i 
to vertex j: 



Wq Wi 
W2 W3 



0|10 



0|00 





llll 



1|01 



Figure 1. Gi for binary CAs with local neigh- 
borhoods of size 2. Each vertex and edge is la- 
belled, the latter by the number to the left of the 
slash. The weight of each edge is the squared 
norm of the amplitude vector of the local config- 
uration indicated to the right of the slash. 



where Wa '■— and a is the base 10 (say) representation of the local configuration 

bit string. The ij^^ entry in A^ is the sum of the weights of the paths with n edges from 
vertex i to vertex j, so it is useful to define the generating function 



n>0 



The i^^ diagonal entry in Ai (t) is the sum of the weights of the paths of length n beginning 
and ending at vertex i, times t"^, for all n > 0; hence TrAi(t) is the corresponding sum 
over all closed paths in Gi. It is straightforward to show [28] that if Z{t) :— det(/ — tA) 
then 

TrAit) ^ (4.1) 
Evaluating the righthand side of (4.1) for Ai we find 



TrAiit) 



{wq + W3 + 2wiW2t)t — 2wQW^t^ 
1 — {wq + 103 + WiW2t)t + WoWsf^ 



(4.2) 



We can read off the weights of the cycles in Gi directly from the positive terms in the 
numerator and the negative terms in the denominator of (4.2): wq, w^, and wiW2- (That 
these are the right terms becomes clear upon expanding (4.2) in powers of t: 



TrAi(t) = {wq + lua)^ + {2wiW2 - 2wqw^ +K + '^wqw^ + wfjt^ + ©(t'^), 

where the cancelling coefficients of have been included to indicate the function of the 
other terms in (4.2).) This is more machinery than we need to find the cycles in Gi, of 
course, but it will become useful in more complicated situations. 



* These are sometimes called 'one-way' automata [15] since the local neighborhood of x extends only 
to one side, although Hillman has pointed out that this terminology is somewhat misleading [34] . 
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Figure 2. G2 for binary CAs with local neighborhoods of size 2. The edges in the subgraph 

isomorphic to Gi are grey and unlabelled; each other edge is labelled by the pair of numbers to 
the left of the slash and weighted by the inner product of the pair of amplitude vectors of the 
two local configurations to the right of the slash. 

For either a periodic or an infinite QCA, each of these cycles must have weight 1 — this 
is condition (z) in the Unitarity Theorems. Thus 

1 = wq (4.3a) 
1 = W3 (4.36) 
1 = wiW2- (4.3c) 

That is, both |00)) and |11)) must have norm 1, while ((10|10)) = ((01|01))-^ 

Figure 2 shows G2 with only the edges in M labelled and weighted; the subgraph iso- 
morphic to Gi is indicated with grey unlabelled edges. Since the orthogonality constraints 
are determined by the weights of paths in M, to simplify the transfer matrix we may set 
the weights of the edges from vertex (0, 0) to itself and from vertex (1, 1) to itself to and 
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the weights of the edges between these two vertices to 1. Then 

(4.4) 



^ woi wol 1 

Wo2 W03 Wi2 Wi3 

Wo2 Wi2 W03 Wi3 

\ 1 W23 W23 



where Wa/3 '■= {{ct\f^)) and a, P are again base 10 (say) representations of the local configu- 
ration bit strings. Now, 

Z2{t) := det{I - tA2) 

= 1 - 2wo3t + (-1 - 2t(;oityo2 + - - 2wi3W23)t'^ 
+ 2{wo3 + W01W02W03 - W01W02W12 - W01W13 

- W02W23 + W03W13W23 - Wi2Wi3W23)t^ 
+ {-'^03 + ^12 + 2woiWo3Wi3 - 2^01^12^13 + 2W02W03W23 - 2wo2Wi2W23)t^ , 

which we may use in (4.1) to conclude that the weights of the parts of the cycles in G2 
which lie in M are wqs, W01W02, W13W23, W01W13 and ^02^23- 

The first two of these are weights of cycles entirely in M, while the last four are 
weights of acyclic paths in M terminating at vertices of Gi. Condition (iii) in the Unitarity 
Theorems forces 

= W01W02 

= W^3W23 
= 1001^^13 
= W02W23- 

Recall that Wa/3 '■= {{a\(3)), so that Wap ~ means that \a)) and \(3)) are orthogonal. Thus 
the only solutions to (4.5) satisfy exactly one of the following sets of relations: 

|00)) ± |01)) A 1 10)) ± 111)) (4.6a) 
|00)) ± |10)) A |01)) ± 111)), (4.66) 

where A is the boolean relation 'and'. These two sets of relations transform into one 
another under the action of P; hence we need only consider one, say (4.6a). 

To obtain a periodic QCA, condition (ii) of Theorem 3.9 must also be satisfied, i.e., 
for the cycles in M: 

= W03 (4.7a) 
= wl^. (4.76) 

Since the amplitude vectors for a binary CA lie in and by (4.3) are nonzero, no more 
than two can be mutually orthogonal. This means that the constraints (4.7) restrict (4.6a) 
to the single frame: 

100)) II |10)) ± |01)) II 111)), (4.8) 
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where by frame we mean a pair Sq, Si of sets of vectors in such that each vector in 
is orthogonal to each vector in Si . We will denote the family of local rules satisfying the 
relations (4.8) subject to the normalization constraints (4.3) by f2^i: the subscript 2 is k; 
the 1 indicates that local configurations i-k/2 ■ ■ - i-ih ■ ■ ■ik/2 with the same state ii have 
parallel amplitude vectors. An explicit parameterization of this local rule gives the rule 
table for /2,l(^|A): 



/2,1 





1 


00 


e*" cos 6 


e'^^i sin 6 


01 


e»</'ipe-*^zsin6' 




10 


gi<^2p-leia cos6' 




11 


e~*^z sin 9 


e~*" cos 9 



where A labels the rows and i the columns. The four entries in rows 00 and 11 form an 
arbitrary SU{2) matrix: a,P,9 G [0, 27r); the remaining degree of freedom is an overall 
phase which has been divided out as it has no effect on probabilities. |10)) is parallel to 
|00)), differing by an arbitrary factor p~^e^'^^ e C. |01)) is parallel to |11)) with length 
p and phase angle (f)i e [0,27r). The only other possible periodic local rule, satisfying 
(4.66), is the the parity transform: /2,-i = -P/2,1, obtained by interchanging the middle 
two rows in (4.9). The state transposition r leads to no new rules: Tin is implemented by 
a ^ it/2 — (3, /9 — s> 7r/2 — a, 9 ^ tt/2 — 9, 0i 02 and p —>■ p~^, while Tout is implemented 
just hy a ^ n/2 + P n/2 + a and 9 7r/2 - 9. 

To obtain an infinite QCA, by Corollary 3.5 and Theorem 3.6, a local rule / must 
have a nonempty deterministic sector consisting of local configurations corresponding to 
the edges in a collection of cycles in Gi. The simplest possibility is that is a quiescent 
state, so that |00)) = (1, 0). Conditions (ii) and (iv) of Theorem 3.10 impose no additional 
constraints beyond (4.3) and (4.5) if 00 is the only local configuration in Df, using (4.6a) 
again we find the local rule /2,i;oo (where the subscript 00 is the deterministic sector), 
which may be parameterized as: 



/2,1;00 





1 


00 


1 





01 







10 




p~^e^^i sin9 


11 




gi<^3e-ia cos 6* 



(4.10) 



since |01)) is orthogonal to |00)), with arbitrary nonzero component pe^'^^ G C, while the 
entries in rows 10 and 11 would form an arbitrary U{2) matrix ^^^^ but for the length 
p~^ of 1 10)): a,/?, 01,03,^ G [0,27r) and the overall phase has been set by the choice 
|00)) = (1,0). Finally, surjectivity constraint (z) of Theorem 3.18 requires in addition that 
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/(0|10) 7^ 0, i.e., cos9 ^ 0; surjectivity constraint (ii) is already satisfied because of the 
orthogonafity relations. Again, the parity transform /2,-i;oo = -P/2,i;00 is also a solution. 

Condition (ii) of Theorem 3.10 will impose additional constraints only if the vertex 
labelled 1 in Figure 1 is also in Di. This occurs if {00, 11} C Df or {00, 01, 10} C Df. In 
either case wi = 1 = W2, since . . . 01 . . . and . . . 10 . . . are admissible configurations in the 
former and . . . 0101 ... is an admissible configuration in the latter. Furthermore, in the first 
case D2nM contains the cycles with weights wqs and in the second case the one with weight 
Wqi, so condition (iv) of Theorem 3.10 imposes constraint (4.7a) or (4.75), respectively, 
either of which restricts (4.6a) to (4.8). This leaves only the completely deterministic limit 
^ = = a = (^1 = 02, p = 1 of the periodic rule (4.9), which we will denote by f2,i-D (the 
subscript D indicates that the rule is completely deterministic), together with its parity 
transform. Injectivity implies surjectivity in this case [32,33], so Theorem 3.18 imposes no 
additional constraints. 

Starting with being an 'anti-quiescent' state, i.e., |00)) = (0,1), forces 11 G Df 
immediately and leads by exactly the same argument to the limit 9 — tt/2 — (3,(f)i — (f)2 = 
0,p = 1 of (4.9), which is again completely deterministic and is immediately identifiable 
as both Tin/2,i;D and Tout/2, i;!?- 

This discussion could be repeated starting with 1 as a 'quiescent' state and would 
lead to the TinTout transformation of the infinite QCA rules found in the previous three 
paragraphs. 

The last possibility for an infinite QCA is for Df to contain only the edges 01 and 10 
which form a cycle in Gi. By condition {iv) of Theorem 3.10, (4.76) must be satisfied, i.e., 
((01 1 10)) = 0, so again (4.6a) is restricted to (4.8) and we are left with exactly the same 
completely deterministic limits already discussed. 

Note that these four deterministic CAs: /2,i;D, Pf2,i;D, Tiiif2,i;D and PTinf2,i;D, are 
trivial in the sense that the local rule depends on only one of the cell states in the local 
configuration, i.e., /2,1;d(«Ki*2) = Su^. That the only deterministic binary CAs with 
neighborhoods of size 2 or 3 are trivial in this sense is well known [35] (this is also the 
sense in which the only linear QCAs which are allowed by the No-Go Theorem are trivial 
for any size neighborhood [8]); we have just shown that if Df contains more than just 00 
(or 11) then infinite k = 2 QCAs are trivial in the same way. It is most interesting that 
when Df = {00} (or {11}) the local rule is not completely deterministic and, furthermore, 
partitions the amplitude vectors into two independent frames rather than the single frame 
of the periodic local rules. 

4.2. k = 3 

The next simplest binary CAs have local neighborhoods of size 3. Figure 3 shows Gi 
(in this subsection the arguments Q — (0, 1} and k — 3 are suppressed). The cycles may 
be determined as in the previous subsection; then condition (z) of the Unitarity Theorems 
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Figure 3. Gi for binary CAs with local neighborhoods of size 3. 
imposes the constraints: 

l = Wo 1 = W2W5 1 = W3WGW5 / , , X 

(4.11) 

1 = W7 1 = W1W2W4 1 = W1WZW6W4. 

Note that only the first five of these constraints are independent: reading the columns left 
to right, the last equation is implied by the three preceding. 

G2 is quite complicated to draw, but by computing successive powers of A2 (now a 
16 X 16 matrix, so computing the determinant in Z2{t) is tedious, even by computer), we 
find the weights of the parts of the acyclic paths of length n < 4 in M terminating at 
vertices of Gi; they are listed in Appendix A. Condition {in) in the Unitarity Theorems 
requires that each of these weights vanish. The only solutions to these constraints satisfy 
one of the following sets of relations: 

|000)) ± |001)) A |010)) ± |011)) A |100)) ± |101)) A |110)) ± |111)) (4.12a) 

|000)) ± |100)) A |010)) ± |110)) A |001)) ± |101)) A |011)) ± |111)) (4.126) 

|000)) II |100)) ± |010)) II |110)) A |001)) II |101)) ± |011)) II |111)) (4.12c) 

|000)) II |001)) ± |010)) II |011)) A |100)) II |101)) ± |110)) II |111)). (4.12(i) 

That any of these sets of relations implies that the weights of all acyclic paths of any 
length in M terminating at vertices of Gi vanish is a consequence of Theorem 5.1 (which 



23 



Unitarity in nonlinear QCAs 



David A. Meyer 



is stated and proved in the next section). Just as in the k = 2 case, the second and fourth 
sets of relations are merely the parity transforms of the first and third, respectively, so we 
need only consider the possibilities (4.12a) and (4.12c). 

To obtain a periodic QCA, condition (ii) of Theorem 3.9 must also be satisfied, i.e., 
we must consider the cycles in M. There are four cycles of length n < 2, with weights: 
wo7, W02W05 and Of the eight cycles of length n = 3, four are contained 

entirely in M, with weights: wosWq^woq, W12W14W24, wirW27W4r and W35W3QW5Q; the other 
four are acyclic paths terminating at vertices in Gi and so are already included in the list 
in Appendix A. Imposing the constraint that each of these weights vanish implies that the 
relations (4.12a) are restricted to a single frame: 

|000)) II |010)) II |100)) II |110)) ± |001)) II |011)) II |101)) II |111)), (4.13a) 

while the relations (4.12c) are restricted to the frame: 

|000)) II |001)) II |100)) II |101)) ± |010)) II |011)) II |110)) II |111)). (4.13c) 

Note that (4.13c) is invariant under parity reversal (and hence is also the consequence of 
(4.12(i)), while (4.13a) is not (and hence (4.126) leads to its parity transform). These three 
frames are the only possibilities for local rules. To show that the longer cycles in G2 rule 
out none of them, consider one and write the transition matrix A2, as in (4.4), with Os for 
the weights which vanish in the given frame. The resulting matrix is sufficiently sparse 
that Z2{t) — det(/ — tA2) can be computed easily. In each of these frames we find that 
Z2{t) = 1 -t^ -2t^ - t^. Since we may compute TrA2(t) from Z2{t) by (4.1), this means 
that there are no further orthogonality constraints on the weights. 

Thus (4.13a) subject to the normalization constraints (4.11) gives a local rule fs^i 
(for k odd we label the local configuration z_(fc_i)/2 ■ ■ - io ■ ■ ■ i{k-i)/2'-i the subscript 1 again 
indicates that the frame in (4.13a) partitions the amplitude vectors according to ii in 
the local configuration). A more concise description of this local rule than an explicit 
parameterization like (4.9) is that the amplitude vectors for /s^i satisfy: 

1 000) ) and |111)) form an orthonormal basis of 
1010)) = ^2|000)), |100)) = ^4|000)), |110)) = ^elOOO)) 
|001)) = ^i|lll)), |011)) = ^3|111)), |101)) = ^sllll)) 

2;j e C,z e {1, . . . , 6}; |2;22;5| = 1, \z\Z2Z4\ = 1, 12:32:5^71 = 1- 

The amplitude vectors satisfying the constraints (4.13c) subject to the normalization con- 
straints (4.11) may be described similarly, changing the middle two lines of (4.14a) to: 

1001) ) = zilOOO)), |100)) = ^4|ooo)), |ioi)) = ^5|000)) 
|010)) = 22IIII)), |011)) = 23IIII)), |110)) = 26|111)). 

We denote this rule by fs^. Finally, the parity transform /3,-i = -P/3,1 gives the solution 
corresponding to (4.126). Each of these rules may be parameterized by 12 real parameters 
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and an overall phase: 6 for the first pair of orthogonal vectors, 2 for each additional vector 
in the frame, and —1 for each independent normalization constraint. These are all the rule 
families for the periodic case; just as when k = 2 the state transposition r leads to no 
additional rules. 

For an infinite QCA, the possible deterministic sectors are determined by the cycles 
in Gi. Viewing Figure 3 or recalling (4.11), we see that the sets of local configurations 
appearing in cycles are: {000}, {111}, {010,101}, {001,010,001}, {011,101,110} and 
(001, Oil, 100, 110}; any union of one or more of these sets is a possibility for Df. 

The simplest possibility is Df = {000}, i.e., is uniquely quiescent. In this case 
|000)) = (1, 0) and conditions (ii) and (iv) of Theorem 3.10 impose no additional constraints 
beyond (4.11) and (4.12). (4.12a) defines the rule family denoted /3,i;ooo with amplitude 
vectors partitioned into four frames: 

|000)) = (1,0); |001)) = (0,^i) 

laO)) and \al)) are orthogonal for a G {01, 10, 11} (4.15a) 
the norms of the amplitude vectors satisfy (4.11), 

a 16 real parameter family of local rules. Surjectivity constraint (i) of Theorem 3.18 
rules out only the codimension 1 submanifolds defined by /(0|010) = 0, /(0|100) = 
and /(0|110) = 0, while surjectivity constraint (m) is again satisfied as a consequence of 
the orthogonality relations. Similarly, (4.12c) would define a rule family with amplitude 
vectors satisfying: 

|000)) = (1,0); |001)) = (zi,0); |010)) = (0,^2); |011)) = (0,^3) 

|10z)) and are orthogonal for i,j G {0, 1} (4.15c) 

the norms of the amplitude vectors satisfy (4.11), 

except that |010)) = (0,^2) conflicts with the surjectivity constraint that /(0|010) = 0, 
ruling out this possibility. Thus /a,!; 000 a-^d its distinct parity transform are the only 
allowed rule families when Df — {000}. Both of these infinite QCAs with uniquely 
quiescent have distinct state transposition transforms under Tin Tout with 1 as a 'quiescent' 
state. 



The next simplest possibility is that Df — {000, 111}. Then by condition {ii) of 
Theorem 3.10, the two acylic paths in Gi terminating at 000 and 111 must each have 
weight 1: 

1 = W1W3 

(4.16) 

1 = W4W6. 

With (4.11), the second of these constraints is implied by the first. By condition [iv) of 
Theorem 3.10, the cycles in D2 fl M must have vanishing weights. Hence = wq^. This 
constraint restricts (4.12a) to the three frame set of relations: 

|000)) II jllO)) ± jOOl)) II jlll)) A jOlO)) ± jOll)) A jlOO)) ± jlOl)) (4.17a) 
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and restricts (4.12c) to (4.13c). We define the rule /3,i;ooo,iii by the set of amphtude 
vectors satisfying (4.17a) subject to the normahzation constraints (4.11) and (4.16) and 
with Df = {000, 111}, described by: 

|000)) = (1,0); |110)) = (^3,0), |001)) = (0,^i); |111)) = (0,1) 

laO)) and \al)) are orthogonal for a e {01, 10} 

the norms of the amplitude vectors satisfy (4.11) and (4.16), 

a 12 real parameter family of local rules. Surjectivity condition (ii) of Theorem 3.18 
restricts these parameter values by removing the codimension 1 submanifolds defined by 
/(0|010) = 0, /(0|100) = 0, /(1|011) = and /(1|101) = 0. Surjectivity condition (ii) is 
ensured by the orthogonality of jaO)) and \al)) for all a e Q^. /3,i;ooo,iii has a distinct 
parity transform but is invariant under TinTout- Applying either of the state transposition 
transforms alone produces a distinct rule family in which is 'anti-quiescent', as is 1. 

The rule which would be defined by (4.14c) with |000)) = (1,0) and |111)) = (0,1) 
is ruled out by surjectivity condition (i) of Theorem 3.18: |010)) = (0,2:2) contradicts 
/(0|010) ^ 0, for example. Similarly, |000)) = (0, 1) and |111)) = (1, 0) fails to be surjective: 
in this case |010)) = (22, 0) contradicts /(1|010) ^ 0. 

We may continue to increase the size of Df and find further infinite QCA rules; the 
procedure is clear. When the deterministic sector gets only a little larger, there will be 
only completely deterministic rules, just as in the k = 2 case of the previous subsection. 

5. Discussion 

The results of the previous section demonstrate that the Unitarity Theorems provide an 
effective procedure for finding one dimensional binary QCAs, both periodic and infinite. 
Although it is increasingly difficult to find the most general unitary solutions for large local 
neighborhood size k — the procedure is not very efficient — the results for k — 2 and k — 3 
suggest a pattern for some particular solutions. Specifically, the sets of relations (4.6) and 
(4.12) generalize to larger values of k. If we generalize the notion of 'frame' to q- frame: a 
collection 6*0, . . . , Sq-i of sets of vectors in such that each vector in Si is orthogonal to 
each vector in Sj for i ^ j, we can state the following: 

Theorem 5.1. Let < j < k e Z. For each 7 e deRne a q-frame 

c(7) _ /c(7) 0(7) 1 

— I'-'O ) • • • ) '-'q-lJ 

by 

Then condition (Hi ) of the Unitarity Theorems — the weight of any acychc path in M(Q, k) 
terminating at vertices in Gi{Q, k) vanishes — holds when the amplitude vectors are parti- 
tioned into the q^ frames {S^'^'*}, or into their parity transforms. 
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Proof. Consider any acyclic path in M{Q,k) starting from a vertex {a, a) e Gi{Q,k) C 
G2{Q, k), where a e Q^~^ . The {k — jY^ edge in the path is necessarily labelled 

{ii - ■ ■ iji'jj^i ■ ■ ■'i'kiH ■ ■ ■ '^j'^'j+i ■ ■ ■ '^'k)i 

where zi . . . ij are the rightmost j states in a and ^ since the first edge of the 
path would not lie in M{Q, k) otherwise. Let 7 = h . . . z^-; then \ii ■ ■ ■ 'ij'ij+i ■ ■ ■ i'D) ^ 'S'|,T^ 
and \ii . . .iji'j^i . . .i'l^)) e S'|7'* . Since all the amplitude vectors in Sp? are orthogonal to 
those in Sp'^ for i" ^ i', 

Hh ■ ■ ■ SS'+i • • • '^fc Ki • • • S'^i+i ■ ■ ■ ^/c)) — 

and the weight of the path vanishes. To show that the parity transformed set of frames 
P{S^~'^} also enforces condition (iii) of the Unitarity Theorems, make the analogous ar- 
gument using the terminal vertex of the acyclic path rather than its initial vertex. ■ 

The special case q = 2 and = 3 of Theorem 5.1 shows that each of the sets of 
relations (4.12) which were found by considering only acyclic paths of length n < 4 implies 
that the weights of all acyclic paths of any length in M({0, 1},3) terminating at vertices 
of Gi({0, 1}, 3) vanish; (4.12a) and (4.126) are parity dual sets of frames for j = 2, while 
(4.12c) and (4.12(i) are the parity dual sets of frames for j = 1. 

That the converse of Theorem 5.1 is false is demonstrated by the existence of the 
nontrivial k = A reversible deterministic CA found by Patt [36] , with local rule: 

1 - Sii2 if h =14 = 0, is = 1; 



f('\'^'^'^'^) = I Su, otherwise. 



Reversible CAs are a fortiori unitary and this local rule partitions the k = 4 amplitude 
vectors inconsistently with each of the sets of frames described in Theorem 5.1. 

Consideration of this example leads to the observation that any reversible deterministic 
CA can be 'quantized': The local rule of such a CA partitions the amplitude vectors 
\ii . . .ik)) into a single q-fraime according to the unique z e Q for which f{i\ii . . .ik) is 
nonzero. Any rigid rotation of preserves this g-frame, and hence unitarity, but gives, 
generically, nonzero transition amplitudes for all the f{i\ii...ik)- The resulting global 
evolution is unitarily inequivalent to the original reversible deterministic evolution. 

Although the local rules for the periodic QCAs found in Section 4 also partition the 
amplitude vectors into a single frame (see (4.8) and (4.13)), they have additional degrees 
of freedom associated with the lengths of the amplitude vectors: 1 when k = 2 and 3 
when k = 3; this should be contrasted with deterministic local rules for which all the 
amplitude vectors have length 1. Despite being asymptotically deterministic, the infinite 
QCAs with local rules found in Section 4 are even further from the deterministic situation; 
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their amplitude vectors lie in more than a single frame: as many as 2 for some of the 
QCAs with uniquely quiescent (see (4.10) and (4.15a)). 

The multidimensionality of the local rule spaces for even the small neighborhood 
QCAs we have considered suggests that binary QCAs may have a wide range of quantum 
behaviors/computational power. Whether any are computationally universal remains to 
be discovered. There is a long standing conjecture that computational power will be 
maximal at critical points of a physical theory [37] . Since the rule spaces here are smoothly 
parameterized this is a more natural arena in which to investigate this conjecture than is 
the deterministic case. 

Consideration of QCAs as physical models, possibly with critical points, raises the 
question of the continuum limits of these models. The simplest nontrivial* one dimen- 
sional linear binary QCAs have the 1 + 1 dimensional Dirac equation as their continuum 
limit [9]. From the perspective of fundamental physics, it would be most interesting to 
determine the continuum limits of the simple nonlinear models we have found here and to 
extend them to higher dimensions. That the reversible deterministic billiard ball model is 
computational universal [38] suggests that higher dimensional QCAs might also be easier 
to prove computationally powerful. It should be noted, however, that there can be no ana- 
logue of the Unitarity Theorems in higher dimensions since reversibility of deterministic 
CAs is undecidable in two dimensions [39]; the best we can expect is, as in Theorem 5.1, 
to find particular sets of local rules which ensure unitarity. 

Despite the existence of computation universal deterministic CAs, probably their most 
important applications are simulations of physical systems [40]. Similarly, it seems likely 
that QCAs will prove optimally suited not to universal computation but for the simulation 
of specific quantum mechanical systems and the solution of particular classes of problems. 
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Appendix A 

The weights of the acychc paths of length n < 4 in M({0, 1}, 3) terminating at vertices of 
G'i({0,l},3) are: 



n = 3: 



n = 4: 



WoiWo2Wq4 

W02WQ4W45 

WoiWo3Wo4Woe 
Wo-i_Wo3WirW26 

W03W04W06W45 

W03WirW26W45 

W04W06W23W47 
W04W23W24W56 

Wo4Wo6W4rW6r 
Wo4W24W56We7 



W01W02W15 
W02W15W45 

Wi5W23W4e 
W15W46W67 

W01W03W06W15 
W01W12W26W35 
W03W06W15W45 
W12W26W35W45 
W06W15W23W47 
W15W23W24W56 
Wo6Wi5W47Wer 

W15W24W56W67 



W01W13W26 
W13W2QW45 
W23W26W57 
W26W57W67 

W01W04W12W24 
W01W03W17W37 

W04W12W24W45 
W03W17W37W45 
W17W23W26W47 
W23W26W35W56 

Wi7W2QW47We,7 

W26W35W56WQ7 



WQ1W13W37 
W13W37W45 
W23W37W57 
W37W57WQ7 

W01W12W15W24 

10011^12^35 W^37 
W12W15W24W45 
W12W35W37W45 
W17W23W37W47 

W17W37W47W67 
W35W37W5eWe7- 



These are determined by (^2)ii (^Dy; respectively, where i,j e {(00, 00), (01, 01), 
(10, 10), (11, 11)}, since this is the set of labels for the vertices in G2 which lie in the 
subgraph isomorphic to Gi. 
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